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Abstract 

We present the calculations of the complete next-to-leading order (NLO) QCD corrections to 
the total cross section, invariant mass distribution and the forward-backward asymmetry (Afb) of 
top quark pair production mediated by W' boson. Our results show that in the best fit point in 
the parameter space allowed by data at the Tevatron, the NLO corrections change the new physics 
contributions to the total cross section slightly, but increase the Afb in the large invariant mass 
region by about 9%. Moreover, we evaluate the total cross section and charge asymmetry (Ac) 
of top pair production at the LHC, and find that both total cross section and Ac can be used to 
distinguish NP from SM with the integrated luminosity increasing. 

PACS numbers: 14.65.Ha, 12.38.Bx, 12.60.-i 



*Electronic address: 



csli@pku.edu.cn 



I. INTRODUCTION 



The top quark is the heaviest particle discovered so far, with a mass close to the elec- 
troweak symmetry breaking scale. Thus it is a wonderful probe for the electroweak breaking 
mechanism and new physics (NP) beyond the standard model (SM) through its productions 
and decays at colliders. The forward-backward asymmetry (Afb) of the top quark pair pro- 
duction is one of the interesting observables in the top quark sector. Within the SM, Afb is 
absent at the tree level in QCD due to charge symmetry, and occur at QCD next-to-leading 

n the last few years, 
lo|. Recently, the CDF 



order (NLO) with the prediction Afb ~ 6% in the tt rest frame [T 
D0 and CDF Collaborations have measured Afb at the Tevatron 7- 
Collaborations annouced that, for the invariant mass of the top quark pair m t i > 450 GeV, 
the measured asymmetry in the tt rest frame is A FB = 0.475 ± 0.1140, which differs by 
3.4cr from the SM predictions Afb — 0.088 ±0.013. This deviation has stimulated a number 
of theoretical papers on NP models, such as new gauge bosons, axi gluons jll - 8 1 1 . 



Recent studies are concerned with the problem of top asymmetry by a flavor-changing 



13 



82j , which can be described by the 



interaction mediated by a charged vector boson, W 
following effective Lagrangian 13): 

Cnp = -gfW'-d-f(f L P L + f R P R )t + h.c, (1) 

where Pr,l — (1 ± 7 5 )/2 are the chirality projection operators, are the chiral couplings 
of the W boson with fermions, satisfying f 2 + f R = 1, and g' is the coupling constant. 
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. It is 



The study of this model at the leading order (LO) has been explored in Refs. 
shown that, for suitable parameters, this model can explain the Afb observed at the Tevatron 
within l-1.5o" of the data. It is well known that the LO cross sections for process at hardron 
colliders suffer form large uncertainties due to the arbitrary choice of the renormalization 
and factorization scales, thus it is necessary to include higher order corrections to make 
a reliable theoretical prediction. Besides, at the NLO level, virtual corrections, real gluon 
emission and massless (anti) quark emission can lead to a sizeable difference between the 
differential top and anti-top production process , which will also contribute to A F b- 
Therefore it is necessary to perform complete calculations of NLO contributions in the W 
model. 

There is a similar work in the Z' model 



84j . where the NLO QCD corrections up to 
0(a 2 g' 2 ) are taken into account. In this work, we calculate both 0(a 2 s g' 2 ) and 0(a s g' A ) 



NP contributions, and the latter term is not definitely smaller than the former so that it 
should not be neglected. Based on the above calculation, we fit the data at the Tevatron, 
including total cross section, the invariant mass distribution and the Ape, and find the 
allowed parameter space. Moreover, we study the top quark pair production at the Large 
Hadron Collider (LHC) induced by a W boson at the NLO QCD level. Since the gluon 
fusion channel dominates in the tt production process at the LHC, it is difficult to probe NP 
effects on Afb from early LHC results. However, LHC will be able to detect the potential 
NP effect on the Charge Asymmetry (Ac) when the integrated luminosity increases in future. 

The arrangement of this paper is as follows. In section [III we show the LO results of 
top quark pair production. In section IHH we present the details of the NLO calculations, 
including the virtual and real corrections. In section [IV] we show the numerical results. 
Conclusion is given in section [V] 



II. LEADING ORDER RESULTS 

Up to NLO, the tt production amplitudes, including NP contributions, can be written as 

M = a sJ S M+9 /np+QWsm + a s9 J N p • ( 2 ) 
And the ti amplitude squared is 

\M*\ 2 = a2/J£/$' + labile {f&f%?) + f^ffi 

+ 2afae (/J£/SH + 2a s V 2 [Tie (&^°*) + ™* (/np /sm°*)] 

+ 2a s g«Tle(f$f$ *). (3) 

In the first line is the LO result of SM and NP. The first term in the second line is the SM 
NLO result and the second term is the interference contribution. The NP NLO result is 
given in the third line. 

The LO Feynman diagrams for the subprocess d{pi)d{j>2) — > t(p3)i(p4) induced by the 
NP and the SM QCD interactions are shown in Fig. [U and the LO partonic cross section 
can be written as 

-LO -LO | -LO , -LO ( A \ 
° = ^SM + ^INT + °NPS I 4 ) 

where subscripts SM, INT and NPS denote the SM channel contributions, the interference 
between SM and NP channels, and NP channel contributions, respectively. The LO partonic 
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FIG. 1: LO Feynman diagrams for dd — > tt. 



differential cross sections are given by 



= -^^[Qmt-AmUt + u)+t 2 + u 2 ], (5) 



d cos 9 9s s 



ddjNT 2/3 a s g 2 / 2 . /-2\ ro 6 , 4 / c 2 . \ 

1 x = 77T — 9 7 — 9 r Jr + Jr 3m? + m. 6m w , — 3t — u) 

dcostf 9sm^s(m^,-t) VJii * v w ; 



+m 2 (t 2 - 2m^, (i + 3m)) + 2m^,« 2 ] , (6) 



,'4 



^^ <,(4-y {(/U/ - )K - 2m ' i 

(4m^, + 4772^,8 + t 2 ) — 8m ( 2 m^ ( tt + 4m^,M 2 ] 
+2/ 2 / 2 [m t 8 - 2m^t + m\ (4m 2 w ,s + t 2 ) - 8m 2 w ,m 4 w , S + Am^s 2 ] }, (7) 

where the Mandelstam variables s, t and u are defined as follows: 

s = (pi+ P2) 2 , t = (p 1 - p 3 ) 2 , u = (p!- p A f. (8) 

The relations between them are 



m 



f = s(l-/9cos0)/2, m 2 -w = s(l + /3cos#)/2, (9) 



where j3 = y/l — Am 2 /s, and 6 is the polar angle of the outgoing top quark in the tt 
rest frame. The colors and spins of the incoming (outgoing) particles have been averaged 
(summed) over. Integrating over cos#, we obtain the LO result of dd — > tt partonic cross 
section. 
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The LO total cross section at the hadron collider is obtained by convoluting the partonic 
cross section with the parton distribution functions (PDF) Gi/A(B) for the initial hadrons A 
and B: 

o" LO = ^2 / dXa I dx b G a / A (x a ,fi f )G b/B {x b ,fi f )a LO , (10) 
where r = 4m^/s. 



III. NEXT-TO-LEADING ORDER QCD CORRECTIONS 

The NLO corrections to the top quark pair production consist of the virtual corrections, 
generated by loop diagrams of colored particles, and the real corrections with the radiation 
of a real gluon or a massless (anti) quark. For the real corrections, we used the two cutoff 
phase space slicing method to subtract the infrared (IR) divergences js5]. 



A. Virtual corrections 



The virtual corrections for the top quark pair production include the box diagrams, 
triangle diagrams, and self-energy diagrams induced by SM QCD and NP interactions as 
shown in Fig. [2] and Fig. [31 respectively. The renormalized virtual amplitudes are given as 
follows 

Msm = T-O e C F (—) M L S ° + [bZ\ + bZ\ + 28g s ) + M f s ™, (11) 

•Mnp = ^C e C F (—) M$ + [5ZI + 5Z\) A^o + M ^ (12) 
47T V e uv/ 

where C e = (47r) e r( - 1 1 _ g . ) . M.(^ and Al^p are ultraviolet (UV) finite terms for SM and NP 
processes. All the UV divergences in the loop diagrams are canceled by counterterms bZ\ for 
the wave functions of the external fields in on-shell scheme, and _5g s for the strong coupling 
constant in the MS scheme modified to decouple the top quark 86], 



6Z% = -^C e ( ) , (13) 



3vr yeuv em 
a s ( 1 2 /i. 



bZ\ = --LC t + — + 4 + 31n^ , (14) 
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FIG. 2: One-loop virtual Feynman diagrams for dd — > tt induced by SM QCD interactions. 
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FIG. 3: One-loop virtual Feynman diagrams for dd — > tt induced by NP interactions. 



where n/ = 5 and fi r is the renormalization scale. The renormalized amplitudes -MgU an d 
A^np are UV finite, but still contains IR divergences. The virtual corrections for subprocess 
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qq — > tt can be expressed as: 



da 



virt 



J ~ virt , J ~ virt , j * virt 
aa SM ' aa NPS ~r ao INT 



2s 



+ 2^e (M%£M™T + M r ^M L s °*) }. 



(16) 



We have calculated the SM contribution, and find the result agrees with that in the Ref. [87] 
The one-loop correction for the cross section induced by NP interactions, with IR singular- 
ities separated from finite terms, is given by 



where 



and 



da 



virt 



+ 



2^ 



a 
a 



(^2)lNT (^i)lNT 
o I" 



(A 



r -ip, 

2JNPS 



em 



ao INT 



+ 



(At 



NPS 



"IP 



em 



dalg s + da mrtJir 



2 J INT 



(A 

(^i)lNT 



-2C F 
4 



rein- 



+ 2 In — - + 9 In ^ + In ^ 



2/9 



In 



//- 

P + l 

1-/3 



20 



(17) 
(18) 

(19) 



(20) 



(^QnPS 

(^i)nps 



-2C F , 



2C f 



2 In 



-ti 



ln4-- 

m t 2 2 



(21) 
(22) 



with t\ — t — mf, U\ = u — ml . The IR divergent terms are proportional to the LO partonic 
crosss section <5"jnt &nd <3"m§q- u v ^ r ^^f^ n jg 

the finite terms of the virtual cross section. 



'nps- 



B. Real corrections 



At the NLO level the real corrections consist of the radiations of an additional gluon or 
massless (anti)quark in the final states as shown in FigHJand FigJSJ 

The phase space integration for the real gluon emission will produce soft and collinear 
singularities, which can be isolated by slicing the phase space into different re gion s using 



suitable cutoffs. In this paper, we use the two cutoff phase space slicing method 



85], which 



FIG. 4: Feynman diagrams for a gluon emission induced by SM QCD interactions. The diagrams 
for a (anti)quark emission can be obtained by crossing the initial-state (anti)quark with the final- 
state gluon. 
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FIG. 5: Feynman diagrams for a gluon emission induced by NP interactions. The diagrams for a 
(anti)quark emission can be obtained by crossing the initial-state (anti)quark with the final-state 
gluon. 

introduces two arbitrary small cutoff parameters, i.e. soft cutoff parameters 5 S and collinear 
parameters 8 C , to decompose the three-body phase space into three regions. 

First, the phase space is separated into two regions by the soft cutoff parameters 5 S , 
according to whether the energy of the emitted gluon is soft, i.e. E 5 ^ 5 s ^/svi/2 : or hard, 
i.e. E 5 > 5 SA /si2/2. Then the collinear cutoff parameters 5 C is introduced to divide the 
hard gluon phase space into two regions, according to whether the Mandelstam variables 
t i5 = (pi — p 5 ) 2 (i=l,2) satisfy the collinear condition \t i5 \ < 5 c Si2 or not. Thus we have 

da Real = da s + da HC + da HV . (23) 

The hard non-collinear term da HC can be written as, 

d^?= J_y |^ 3 | 2 rfr 3 |^ (24) 



which can be evaluated numerically using standard Monte-Carlo techniques 



In the 



following sections, we discuss the parts containing the soft and hard collinear sigularities. 



S 



In the limit that the energy of the emitted gluon becomes small, i.e. E$ ^ £ SA /si2/2, the 
three-body cross section da s can be factorized as 



da- 



Ms 

V2k 



^INT /?NPS 
LO ^ij I ^iLO 
INTENT "° 







XPSqNPS 



-Pi ■ Pj 



(Pi ■ Ps) (Pj ■ Pb) 



(25) 



where C' e = ^(1-2% ("sir) • ^ e c °l° r charge factors Cij are 



/7INT 
L '12 



/i INT 
u 34 



C F /2, 



/?INT _ ^> INT _ ^2 
u 33 ~~ L '44 ~~ ^A^/?, 



/7INT _ /7INT _ _^ ^2 



13 

INT 



c. 



24 

INT 
22 



/7INT _ /7INT _ n IO 
u 14 ~~ u 23 ~~ I 



(26) 



and 



£NPS _ /7NPS 



'13 



/?NPS _ ^7 

u n — u 



U 24 — OO/IOf, 



nps 

22 



0. 



/7NPS _ /^NPS 
u 33 ~~ u 44 



C 



SCaCf, 
nps /7nps /7nps /7nps 



12 



31 



1-1 



a 



23 



(27) 



Here, Cq NT — CaCf and C^ PS = 9 are the color factors o 



'AS^F 



' LO diagrams in Fig. [TJ 



The integration over the soft phase space is given by 



85j: 



1 / 4 \ " e 

dS=-[ — / dE 5 E}~ 2e sin 1 ' 26 9 1 d9 1 x sin" 26 9 2 d9 2 . 

k \S12J Jo 



(28) 



We define 



Then we have 



T. 

-£--7 



(Pi ■ P5) (Pj ■ P5) ' 



(29) 
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Z12 

X14 



Z22 
X44 

2 

s 
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s 



' ' 1 q-f in 
"rr,2 I +X 33 > 

m f <tr 



em /3 1-/3 



4- + — f-21n<y. 



'34 5 



X 



21 



z 



23 



1 
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TR 
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1 

em 
1 
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1 q-fin 

' x 12 ) 

2 In — + ]xiAr + 2la5 i 



7U 



2 In 



^I + ln^ + 21n5< 

s mf 



1 q-fin 
' x 13 ) 



1 q-fin 



(30) 
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where all the IR sigularities in have been extracted out and for briefness, the finite terms 
X(p are not shown here. Now, the Eq. (j2"3"|) can be rewritten as 



da s = ^C, 



+ 



=IR 



em 



+ « 



INT 



da LO 

au INT J 



{A'2 )nps (Af )nps 



IR 



em 



+ (^o)nps 



tt0 NPS> 



(31) 



in which 



(Af)lNT 
{A x )lNT 



2CV, 
1 



161n-^ + 21n-^ + 91n-^ + ln^ 



and 



(A 2 ) NPS 



(Af 



NPS 



+ 161n5 s - 

2C F , 
-2C F 



K 

(1 + /3 2 ) 
2/3 



//2 



mr 



In 



1-/3 



(32) 



21n-^ + ln-^ + 21n5 s -l 



(33) 



In the hard collinear region, ^5 > £ sv /si2/2 and |t i5 | < 5 c Si 2 , the emitted hard gluon is 
collinear to one of the incoming partons and the three-body cross section is factorized as 



da 



HC 



da 



LO 



0:.s 



L2^ 



- ] 5 c e [Pdd{z,e)G d/p (xi/z)Gd/ p {x 2 ) 

dz ( 1 — 



+PSd{. z ^) G d/ P {^-i./z)G i / v {x^) + ( Xl ^2)] 



dx\d%2i (34) 



2 \ 2; 

where P^ are the unregulated splitting functions in n = 4^2e dimension for < z < 1, which 
is related to the usual Altarelli-Parisi splitting kernels 
Explicitely, in our case, 



89(] as Pifae) = P^z) 



ePlj(z) 



Pdd(z) 



Pdd(z) 



Pdd(z) 



C 



1 + z 2 



1-z' 
C F (l-z). 



(35) 
(36) 



For massless d(d) emission, we decompose the phase space into two regions, collinear and 
non-collinear, and give the expression for gd — > tid cross section, 

da(qg — > tiq) = a°(ag — > tia)[G a / p (x i)G g /p(x 2 ) + (x 1 <r» x 2 )]dxidx2 

ce=d,d 



+ da 



LO I a s n l 



1 



— ^c e [ p d 9 {z,e)Gg /p {x 1 /z)Gj /p {x2) 



dz (\ 

+ P-[ g (z,e)G g / p (x 1 /z)G d / p (x 2 ) + fa x 2 )} — 



dx\dx2i (37) 
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where 



P dg (z) = P-Az 



and 



dg\ 



o c [ag — > tta) 



W + (l-z) 



2s 12 



K 9 {z) = P' dg (z) = -z(l - z) 



\M 3 \ 2 (ag -)■ tta)dT 3 \ H ^. 



(38) 



(39) 



In order to factorize the collinear singularity into the PDF, we introduce a scale dependent 
PDF in the MS convention |85j], 



G a / p {x,Hf) = G a/p (x) + ^2 (~^) 



2w \ 



As in Ref. [85J], the 0{a s ) collinear contribution is 



1 dz 

— P ,p{z)Gp /p {x/z). 

z 



da 



coll 



where 



+ £ 



a=d,d 



{Gd/ P (xi, Hf)G d -/ p (x 2 , /i/) + G d / P (x 1 , fif)G d / p (x 2 , [if) 



+ Ao c (a — >■ ag) 



Gd/p (x i ,fif)G d / p (x2iHf) 



+ (xi -H- X2)}rfXidX2, 



(40) 



(41) 



and 



with 



(d ->• d#) = (d -> d#) = C F (25 S + 3/2), 



= Al c ]n Sl2 



Ga/p ( % 



a J X 



2 ' 



l-SsSa/3 



dy 



G/3/p(xi,Hf)P a ^(y), 



Finally the NLO correction of dd — > ti process can be written as 



a 



NLO 



(42) 



(43) 



(44) 



dx x dx 2 [G d/p (xi, ^/jGj/pCxa.M/) + (&i ^ £ 2 )] (<7™ rt + a 5 + a HC ) + o coU j 
+ ^ ^ dx x dx 2 [G g / P (xi, fif)G a / p (x 2 , fif) + (xi <H- x 2 )] o c {ga -> tta). (45) 



a=d,d 
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Note that all the IR divergences in the NLO total cross section are proportional to the 
LO cross sections, and we find the following relations 

(^)int + (Af ) rNT = 0, (Al) mT + (Af) mT + ^ Af{a -> ag) = 0, 

a=d,d 

(^)nps + (A§ ) NPS = 0, K) NPS + (Af ) NPS + J] ^i C (« -> ag) = 0. (46) 
Now all the IR divergences are canceled exactly. 

IV. NUMERICAL RESULTS 

In the numerical calculations, we set mw> = 400 GeV, because such a W is readily 
observed at Tevatron with an integrated luminosity of 10 fb _1 , and at the LHC with an 
integrated luminosity of 100 pb _1 [32|. There are two independent parameters in the NP 
Lagrangian. For the convenience of calculations we define the a parameter set (Cy, Ca), 
where Cy = g'(fR + fiJ/2 and Ca = g'(fR — fiJ/2. The mass of top quark is chosen to 



be m t = 172.5 GeV. The CTEQ6L and CTEQ6M PDF sets 90] and the associated a s 
functions are used for LO and NLO calculation, respectively. Both the renormalization and 
factorization scales are fixed to the top quark mass unless specified otherwise. 



A. scale dependence 

In Fig. [6] we show the scale dependence of the LO and NLO total cross sections at the 
Tevatron for three cases: (1) the renormalization scale dependence \i r = fi, [if = m t , (2) the 
factorization scale dependence /i r = m t , [if = //, and (3) total scale dependence 
From Fig. [6j we can see that the NLO corrections significantly reduce the scale dependence 
for all three cases, making the theoretical predictions more reliable. 



B. Tevatron constraints 



A-p-Q of top quark pair productions is defined as 



A 



FB 



Op + (T B 

A™xR + Af£x(l-R) 



12 



f,=o,f„=i 



\i=[L, \x=m t 
\i=\i, \x=m t 




FIG. 6: Scale dependences of the total cross sections at the Tevetron. The black and the red lines 
represent the LO and NLO results, respectively. 

where 



.4 



NP 
FB 



cr 



NP 



a 



NP 



ASM 
^FB 



R 



SM 



a 



SM 



)/(* 
)/(* 



NP 
F 

SM 
F 



+ °B ), 
i ^SM\ 



a. 



_NP 

tot 



SM , NP^ 
tot ' °tot J 



(47) 



are the asymmetries induced by NP and SM, and R is the fraction of NP contribution to 
the total cross section, a-p and <7b denote the total cross sections in the forward(F) and 
backward(B) rapidity regions, respectively. The LO and NLO total cross sections of the 

v and C 2 A , 



interference and NP contributions can be written in terms of Cv and C 2 



INTl 

a Jlo 



:i-i4)i^(c v + ci)] P b, 



Cr 



NPS1 

tt Jlo 



-0.31 V 

[2.06183 W + C\f - (2.51)11^(0* ■ C|)l pb, 



(48) 
(49) 
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FIG. 7: Values of Cy and Ca allowed by Tavetron data at 95% CL: <7 t j=(7.50 ± 0.48) pb , 
A FB (mtt > 450 GeV)=0.475±0.114, and (^ t 7dm< f ) m " e[800 ' 14001 GeV = (0.068 ± 0.036) fb/GeV. 
The blue dot (0.78, 0.78) and brown star (0.74, 0.74) represent the BFPs at LO and NLO level, 
respectively. The allowed parameter region is symmetric with respect to the Ca and Cy axes, so 
we only display the contours where CaCv > 0. 



and 



[°ti J NLO 

.NPS] 
tt J NLO 



[-(1.42)^)(CHCi)] pb, 

t 



2.39^(C 



C A ) 2 



(2.82)^ 4 (C v • Ci 



pb, 



(50) 
(51) 



where the errors are obtained by varying the scale between \x r = \if = m t /2 and fi r = fif = 
2mt. The differences of the cross sections in the forward and backward rapidity region are 
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given by 



(J 



INT 



cr 



NPS 



a 



a 



INT] 
B JlO 

NPS] 
B JlO 



[-(0.26) + ^(C v + CD] pb, (52) 
[0.68;°;i°(C v + C\f + 0.054tS(C v ■ C\)} pb, (53) 



and 



r^.INT JNT] 
L°F a B J NLO 



[« 



NPS 



cr 



NPS] 

B J NLO 



[-(0.40);°;°4(C* + CD] pb, 
[0.94 



0.03/n 2 i p2\2 / n 107 \ -0.000/^2 
+0.04 vW ~r ^A) — l u - iz ' J+0.0021W 



Ci)] pb. 



(54) 
(55) 



For tt invariant mass spectrum, we restrict our attention in the large invariant mass region, 
i.e. m t t € [800, 1400] GeV, where the A FB is the most obvious. The results are presented as 



and 



. dmu 



- d*r 

. dmu 

d*r s 

dmu 



m tf £ [800, 1400] 
LO 

m tf £ [800, 1400] 
LO 



m t t £[800,1400] 
NLO 

m tf £ [800, 1400] 
NLO 



Ro.oi4)l» v + ci)] 



pb 

GeV' 



[0.082i°;™(C v + Ci) 2 - (0.064)±85£(C^ • Ci)] 



Pb 

GeV 



(56) 
(57) 



[-(0.012)i°;^(C v + C 



pb 

GeV' 



[0.117^(C V + Ci) 2 - (0.094)^(C V ' Ci 



pb 

GeV' 



(58) 
(59) 



From the errors in Eqs.fl48l- 159]) we can see that NLO corrections reduce the dependence of 
the cross sections on the renormalization and factorization scales. 

In Fig. [71 we show the allowed region in the (Cy, Ca) plane that is consistent with the 
total cross section a t t, Afb [9] and the spectrum of m t i in the large mass region 91 1, which 
are given by 



a 



EX 



day. 
dm+ 



a EX 
FB 

-1 m tt -£ [800, 1400] GeV 



(7.50 ± 0.48)pb, 

0.475 ± 0.114, for m& > 450GeV, 
(0.068 ± 0.036) fb/GeV. 



(60) 
(61) 

(62) 



EX 



We use Monte Carlo programm MCFM [92] to get the cross section of the gluon fusion 
channel gg — > tt at the NLO QCD level. As for the process of qq — > ti, we have checked our 
value with the results given by MCFM at QCD NLO level, which are well consistent in the 
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range of Monte Carlo integration error. Combining the contributions of these two channels 
we have the SM predictions for the above observables at NLO level 

4 M = 7.00±g;?g pb, (63) 

m tt -e [800, 1400] GeV 

= 0.055™ fb/GeV, (64) 

SM 

where we have considered scale uncertainty in the calculations. For consistency we have 
used the SM QCD predicted values of A FB (m tl > 450GeV) = 0.088 ± 0.013 at NLO level, 
although next-to-next-to-leading logarithmic (NNLL) SM QCD results are available [6]. 

The measurements of /Lfb and invariant mass spectrum da t t/dm t i in the large invariant 
mass region are particularly sensitive to values of Cy and Ca at NLO level. In order to 
generate the desired A F b m large m t t region, NP couplings should be large enough so that 
the positive NPS terms could overcome the negative INT terms. While on the other hand, 
the NLO NPS effect causes the cross section in the last bin of m t i to exceed the la upper 
limit of experimental result and therefore we expect the couplings Cy and Ca to be not too 
large. As a consequence, NP couplings are subject to strong restrictions. 

In Figs.[7J solid and grid regions correspond to NP LO and NLO results at 95% confidence 
level (CL), where we have considered theoretical and experimental uncertainty in a t i and 
dattjm. t i and only consider experimental uncertainty in the Afb calculation. It can be seen 
that NLO corrections manifestly trim down the area of the allowed parameter region. The 
blue dot (0.78,0.78) and brown star (0.74,0.74) represent the best fit points (BFPs) at LO 
and NLO level, where x 2 reaches its minimums of 2.0 and 2.3, respectively. Thus we can 
see that higher order corrections impose stronger constraints on NP couplings and reduce 
the BFPs of C v and C A by 5%. 

Now we discuss the theoretical predictions for the measurements at Tevatron induced by 
NP at the NLO BFP (0.74,0.74), or equivalently, g' = 1.48, and f R = 1, f L = : The LO 
and NLO total cross sections of tt production are 

< p LO = 0.461 pb, 

4,nlo = 0-458 pb, (65) 
and the differential cross section are 

[rf4 P /^-]^ ie[800 ' 1400]GeV = 0.062 fb/GeV, 

[rf4 P /^-]^o [800 ' 140 ° 1GeV = 0.083 fb/GeV. (66) 
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da t i 
dm t t 



Here, the superscript "NP" represents the combination of the "INT" and "NPS" contri- 
butions mentioned above. It can be seen that the NLO corrections have slight effects on 
the total cross section but increase the invariant mass distribution in the large mass region. 
Note that the two parts of the NP corrections, INT and NPS terms, are individually not 
small, but they have opposite sign and cancel each other. 

The Afb containing NP contributions at the NLO BFP are shown in Table [H All the 
theoretical predictions containing NP NLO effects are consistent with experimental results 
within 2cr CL. It is found that the Afb in the large invariant mass region gets an obvious 





SM NLO QCD + NP LO 


SM NLO QCD + NP NLO 


APP 
^FB 

Att 
^FB 

A% (m a < 450 GeV) 
^fb ( m tt > 450 GeV) 


0.130 
0.146 
0.074 
0.250 


0.140 (0.2 a) 
0.156 (0.0 a) 
0.068 (1.2 a) 
0.272 (1.8 a) 



TABLE I: The A FB with £ = 1.48 , f R = 1, f L = and M w > = 400 GeV at the Tevatron, where 
A F g and Ap B are the Ap^ in the lab frame and the tt rest frame, respectively. Here we list the CL 
when including NP effects at NLO level. 

enhancement by about 9%. 

In Fig. [HI we show differential cross section da/dm t t when we consider NP effects at 
the NLO BFP, from which we can see that higher order corrections do not change the 
distribution very much. 



C. LHC predictions 

The process of top quark pair production has been measured at the LHC, and the cross 



section 



93 



94| is: 



c4 TLAS (LHC, VS = 7 TeV) = 180 ± 18 pb, (67) 
ag MS (LHC,VS = 7TeV) = 158 ± 19 pb, (68) 

which are consistent with the SM predictions. Since the LHC is a proton-proton collider, 
which is forward-backward symmetric, the Afb defined at Tevatron can not be directly 
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FIG. 8: Differential cross sections da/dm t i as a function of m t j at the NLO BFPs (±0.74, ±0.74). 
Here " Ex per imental data" is da/dm t i measured with 2.7 fb _1 of integrated luminosity at the 
Tevatron |9l|. "SM NLO QCD" represents the results in the SM QCD at NLO level. "NP LO + 
SM NLO QCD" and "NP NLO + SM NLO QCD" stand for the predictions including NP effects 
at LO and NLO level, respectively. 



applied to the proton-proton collider experiments at the LHC. Ac used by CMS [95|, |96| is 
denned as 



A f 





- 


Vt\ > 0) - cr(|77 t 


-h 


<0) 






m 


>0) + a(\ Vt 


- h 


<0) 



(69) 



where rj t and r\i are pseudo rapidities of top and antitop quark, respectively. Last year at 



CMS, it is measured to be [95[], 



A c = 0.060 ± 0.134(stat.) ± 0.026(syst.) 



whereas the recently updated report shows 96j 



A c = -0.016 ± 0.030(stat.) ± 0.019(syst.) 
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FIG. 9: Results of a combined fit to er t j and the value of Afb allowing for NP at different CL. 
The shadows from dark to light indicate, the experimentally favored region of 95%, 97%, and 
99% probability in the Cy — Ca plane. The black dashed lines and grey solid lines, respectively, 
represent the value of the total ti cross section and the Ac at the LHC with y/S = 7 TeV. The 
black dot and stars represent the SM point and the NP NLO BFPs . 

The discrepancy between these two mesurements is evident. However, given the large ex- 



perimental error, both results are compatible with the SM predictions Ac = 0.013 95]. 
The A c induced by NP interactions at the NP NLO BFPs (±0.74, ±0.74) is 0.069, which is 
about 5 times of SM prediction. This result is very close to the previous observed central 
95], and is also consistent with the the latest data value {96] within 2 a CL. 



value at CMS 



We still need more experimental data with higher precision to seek evidence for a possible 
modification in Ac by NP. 
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In Fig. |9l we show the results of a combined fit to the ti data in the presence of NP at 
different CL. The shadows from dark to light indicate the experimentally preferred region 
of 95%, 97%, and 99% probability in the Cy — Ca plane. The black dot represents the 
SM point (0,0), and the black star represent the NLO BFP (C V ,C A ). At the LHC with 
VS = 7 TeV, the cross section of tt production at NLO QCD level in SM is a t t = 154.5 pb. 
Including NP contributions at the NLO BFP, we have oJgSJ 1 = 175 pb. 

From the shape of contours of a t i and A c , one may easily distinguish NP events from SM 
ones. The location of the brown contours indicate that vector current Cy or axial currrent 
Ca alone cannot improve the quality of fit significantly. The acceptable confidence region 
(global 95% CL) centers around four points in the parameter space where |Cy| equals |Ca|, 
which means experimental data favor purely right-handed or left-handed couplings. 



V. CONCLUSION 



We have investigated NLO QCD effect on total cross section, invariant mass distribution 
and forward-backward asymmetry Afb of top quark pair production mediated by W at the 
Tevatron and LHC. We have taken into account the interference of NP channel with QCD 
channel ( up to 0(a 2 g' 2 )), as well as the interference between NP channels( up to 0(a s g' A )). 
We fit the data at the Tevatron, including total cross section, the invariant mass distribution 
and the Ape, and find the allowed parameter space. We show that due to the cancelation 
between these two parts of contributions, the NLO total cross section exhibits only a slight 
modification compared to the LO result of NP. But the Afb is increased by about 9%. 
Moreover, NP couplings is constrained strongly by the discrepancy between NP and SM in 
the invariant mass distribution. And these constraints become more stringent in the QCD 
NLO. Thus, it is difficult to satisfy simultaneously the constraints from data of Afb and 
da/drritt spectrum in the large invariant mass region at the Tevatron. At the LHC, both 
total cross section and Ac can be used to distinguish NP from SM and therefore the LHC 
may detect these NP effects with the integrated luminosity increasing. 
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